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ABSTRACT 


Binary sequences have had application in communication 
Systems tor many years. Shift registers have been used in 
their generation, because of the ease and economy of their 
Operation. For certain applications, nonlinear feedback 
functions are used by shift registers of span n to generate 
sequences of lengths up to oe 

The sequences of maximum length 2" and their generation 
are the subject of this thesis. [In particular the ways 
Of generating these sequences using nonlinear feedback 
shift registers and their correlation to linear feedback 
shift registers are described. Complexity is the term 
fmmen to the lemgth of the shortest linear feedback shift 
register generating a maximum length as sequence, 

Games and Chan [Ref. i] have given considerable study 
to the subject of complexity. Some of the problems they 
left are discussed further in this paper. It will be shown 
that the complexity of a de Bruijn sequence (S) is the same 
as the complexity of its reverse (r S), complement (S), and 
its reverse complement (r S). Sequences (S) for which r S = 
S are termed RC sequences. It is shown that RC sequences 
Pest £On every Odd n>5. in addition a lower bound will be 


established for the number of RC sequences occurring for 


each odd n. 
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fm oN ROwUs 1 TON 


fw APPLICATION 

de Bruijn sequences have been the subject of recurring 
interest since their discovery over a century ago. Recent 
developments in radar ranging [Ref. 2: Ch. 13], error correct- 
ing codes, secure or limited access code generators [Ref. 3: 
pp. 12-14], and mathematical modeling have made application 
of de Bruijn sequences because of their properties. The 
binary shift registers readily model digital computers with 


electronic states where "l'' indicates the on state and ‘'0Q" 


indicates the off state. 


PeeelEFINITIGN OF de BRUIJN SEQUENCES 
For a positive integer n, a de Bruijn sequence of span 
ets aacomplete baianary cycle of length 2” which is a sequence 


Se= 45.5 s,, °*-, s _} taken in circular order (s, follows 
i Z on i 
S F) such that all possible n-tuples occur exactly once. 
2 
[Ref. 4: p. 120]. Complete binary cycles are often called 


de Bruijn sequences after the Dutch mathematician, N. de 


Bruijn, who in 1946 proved the existence of complete binary 
; mt 
cycles having length 7 ad tat they numbered ae as 


de Bruijn was preceeded by Flye Saint-Marie a half century 
earlier in an obscure publication [Ref. 5]. As an example 


2°13 


for n = 3, Table 1.1 shows both (2 ) de Bruijn sequences 


ll 





of length eight ar Born OL these sequences are examined 
MOmeOiorrm ti@emeach Of the eight 3-tuples occur exactly 
Miee. NOtice that 2 elements of the sequence are repeated 
as signified by the bar (00), which was necessary on each 


sequence to complete the last two 3-tuples. 


TABLE 1.1 
Sete ee COMPOSITION OF de BRUIJN SEQUENCES OF SPAN 3 


0001110100 0001011100 


000 000 
001 001 
Seo 010 
I 101 
110 O11 
101 vie 
010 ee 
100 100 


There are 2” different possible starting points for 
each de Bruijn sequence; however, each cyclic permutation 
1s considered equivalent. It is often a matter of conveni- 
ence to start each sequence in a canonical way with the all 
zero or all one n-tuple. This reduces confusion when working 


with a number of different de Bruijn sequences. 


C. PSEUDO-RANDOMNESS 

Since each binary n-tuple occurs exactly once on the 
cycle, the sequence models a uniform distribution with the 
probability of given n-tuple occurring to be exactly = ee 


The sequence is not entirely random, however, since each 


Ly 





n-tuple (ay a a) has but two possible successors 


(0 a, a, ce* a or l a, a, *** a,_j)- de Bruijn sequences 


ie ik 
Satisfy the randomness properties of equidistribution of 
0's and 1's and run lengths [Ref. 3: p. 10] which would be 


expected from the tossing of a fair coin. Thus, de Bruijn 


sequences have a pseudorandom property. 


tee &-ARY SEQUENCES 

In general it 1s possible to work with k possible states 
for each position of an n-tuple. Flye Saint-Marie in 1894, 
Showed the existence and determined the number of complete 


,a-l ,n-l_y 
cycles for n-tuples of k characters to be [(k-1)!] k : 
[Ref. S]. This report is restricted to the binary case 
fiow= 2) due to its principal application in electronic con- 


munications and computers. 


E. DEFINITIONS 
For S = {s, Eo°** S;} where S. € CORE forme i<i tiie 
following are defined: 
tL. Let @(S) denote the weight of S. 
k 
W(S) = = S5 
i=] 


7. beter to) demote the parity of S$. 


P(S) = gl (mod 2) 





cect lio) demote the length of S. 


1(s) = k where k is the number of positions in S. 


14 





IIT. SEQUENCE GENERATORS AND OPERATORS 


Pe SHIFT REGISTERS 

In practice shift registers are used to generate de Bruijn 
sequences. An n-stage shift register has n memory registers 
(Xx), Xor °%%, x) which shift their contents (0 or 1) to the 
Mext register upon command. For example, in Fig. 2.1, the 
contents of x4 will transfer to Xo, Xo to Xz, etc., at the 
appropriate time with serving as the output. However, 
tme contents of an n-stage shift register would empty in 


n shifts if no input was provided to Xj: 


OUT 





Pee, 2s n-secage Shift Register. 


Ea, PEED ACK SHIFT "REGISTERS 

Feedback networks are added to provide an input to the 
met stage of @he shaft register. Note that in practice 
the output can be taken from any stage of the shift register. 
eect incecontents Of @miy stage is the same as the con- 


Wemto OL any Stage is the s@me as the contents of the nth 


15 





eee wmerely shifted by a certain amount. In Fig. 2.2 the 
feedback function f is added to the n-stage shift register 
to generate non-trivial sequences. The coefficients 


ag, 44, °**, a, are 0 or 1 with ay = a, = le, 
a 


Feedback 
Loop 


OUT 






Fol Oteecee n-Stage Feedback Shift Register. 


1. Linear Feedback Functions 
A feedback function f is classified as linear if 
the function f (x) > Xo > rb eree.) is restricted to be of the 
form £ - a, X; where the addition is modulo 2 addition. 
The Sinisa (mod 2) table is given in Table 2.1 below. 
Note that subtraction is equivalent to addition when opera- 


ting modulo 2. 
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ee 2a.) 


ADDITION (MODULO 2) 





2. Nonlinear Feedback Functions 
If the feedback function E(x, , coe, x) utilizes 
multiplication (mod 2) and addition (mod 2) then the feed- 
back function is nonlinear. (We also say that the linear 
feedback functions are vacuously nonlinear feedback functions. ) 
Modulo 2 multiplication is identical to multiplication in 
the integers with the restriction of the alphabet used 


being the set {0, 1} as described in Table 2.2. 


TAbHeEY 2.2 


MULTIPLICATION (MODULO 2) 





3. Comparison of Linear and Nonlinear Feedback Functions 


Considerable research has been devoted to the study 
of linear feedback functions because of their ease of analy- 


Sis. Nonlinear feedback functions are much more complicated. 


ey 





Hovever, the class of sequences that can be generated by 
nonlinear functions 1S much greater than for linear functions. 

The number of binary maximum length (2"-1) linear 
sequences that can be generated by a n-stage linear feed- 
back shift register (LFSR) is approximately Li [Ref. 7]. 
The reason that the maximum length linear sequences only 
have length 2™-1 is that the zero state is a fixed pormnt 
under the linear function. This poses no problem since the 
maximal length linear sequences are readily made into 
de Bruijn sequences having length a Simply by adding a 
zero to the (n-l)-tuple of zeros. The feedback function 
involved then becomes a nonlinear function. 

By comparison the number of binary nonlinear full 


length sequences a tiat Gan oe generated by a nonlinear 


Mivl 
feedback shift register (NFSR) is exactly ye wa 


Below, 
Table 2.3 compares these numbers for 3<n<8. The number of 
maximum length sequences will be of interest later in 

section III, where the concept of complexity in generating 
these sequences will be developed. Various algorithms are 


listed in a survey by Fredricksen [Ref. 5] for generating 


age Bruijn sequences. 


Gey SUCCESSOR AND PREDECESSOR STATES 
The contents of an n-stage shift register at a specific 
time 1S usually referred to as its "state". As mentioned 


in Section I when discussing the pseudorandomness of de Bruijn 


18 





Pwo, 2 «5 
NUMBER OF MAXIMUM LENGTH SEQUENCES 


n Linear Nonlinear 
3 2 2 

4 4 16 

: 6 2048 

6 6 67,108,864 
7 18 Le ve Tad 
g 30 7 hate” 


sequences each state has two possible successor (conjugate) 
States depending on whether the feedback function generates 
foe )6COr CL" «as the input to register X4- Likewise each 
State has two predecessor (companion) states depending on 
Miaethner the previows output was a "0" or "1". The adjacency 
G@u@@adruple in Fig. 2.3 centers around then (n-1)-tuple Be °c 
That is, each input state has two possible successors depend- 
ing on whether the feedback is 0 or 1, and each output state 
has two possible predecessors depending on whether the pre- 


wh Stage was 0 or l. 


vious output from the n 
In Table 2.4 we show an example for n = 3 of the possible 


successors and predecessors for each state. 


1? 





Fig. 2.3 Adjacency Quadruple. 


TABLE 2.4 
CONJUGATE AND COMPANION STATES (n=3) 


Predecessors State Successors 
000-001 000 000-100 
010-011 001 000-100 
100-101 010 001-101 
tO =e 011 001-101 
000-001 100 010-110 
010-011 101 010-110 
100-101 110 011-111 
110-111 He i 011-111 

a +7 
Output Input 


D. de BRUIJN DIAGRAMS 
The de Bruijn diagram compactly contains all the above 
information. The diagram contains 2" vertices corresponding 


to the 2” states and two directed edges from each state to 


20 





the possible pair of successor vertices. As a result two 
arrows exit every vertex leading to conjugate states, and 
two arrows enter every vertex arriving from companion states. 
For an example consider the de Bruijn diagram Go for n = 3 


im Fig. 2.4. An Eulerian path through G is defined as a 
010 
101 


110 Q11 


dela) 





Patiewhich Visits each edge exactly once. A path which 
Visits every vertex of oe 1s called a Hamiltonian path. A 
Hamiltonian path in G, 1s a de Bruijn sequence of span n. 
[If the edges of G_ are labeled with an n + 1-tuple defined 
by the labels on the predecessor and successor states, then 
the Eulerian paths in = correspond to Hamiltonian paths-- 


and de Bruijn sequences in the graph Go The Hamiltonian 


+] ° 


Maun in Fig. 2.5 is a subgraph of G. in Fig. 2.4 with exactly 


one edge emanating from each vertex. 


E. OPERATORS 


For S = S$; S5 °°" s where S; mel ton l<i<n the 


S 
n-l “n 


following operators are defined. 
i identacy Operator 
The identity operator (e) operating on S denoted 


eS = Sy S5 eee a oe 


-n-1 n 
2. Reverse Operator 


The reverse operator (r) operating on S is denoted 


r S, where r S = s_ 5s 


= eae S5 S,;- The reverse operator 


preserves the weight, parity and length of S. Note also 
that (r)* =e. 
5. Complement Operator 
The complement (dual) operator on S is denoted d § 


am 5, where & = Foo a, Se such that Ss. = s, 91 


(mod 2) for l<i<n. 





000 


ale 100 
010 
101 

O11 110 
hdd 


Fig. 2.5 A Hamiltonian Path through G.. 
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beme 1: If W(S) = kefor S = S$; 5S, °°" S then 
W (S) = n-k. 
Boot: W (Ss) = kK implies k nonzero elements ec S 
amd n-k zero elements _« S. 
Complementing n-k zero elements produces 


n-k nonzero elements ec S. 


W (S) = n-k 
One 
Lemma 2: If P (S) = a where a c« {0,1} and S = {s, Sy °° See, 
then P (S) is a iff n is even. 
n 
Proof: P (S) = 2 S5 (mod 2) =a 
1=1 
ie n 
P (S)e2R Er S; (mod 2) = ¢f S5 ®l=ae@n 
1=1 1=1 
therefore P (S) = P (n-a) 
Orie). 


Also note ane =e, and that the reverse and comple- 
ment operators are commutative, i.e. rd = dr. 


4. . Reverse Complement Operator 


Wac reverse complement operator oneS is denoted 


rS= Ss, S,.1 °°* $2 8, - Weight and parity are effected 
by the reverse complement Seema’: exactly as with the 
complementation operator. 

These operators apply to sequences as well, since 
Sequences are composed of n-tuples. Table 2.6 shows the 


Pmiect Of Opemators on de Bruijn sequences of span 3 and 4 
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where r S, S, and r S are all de Bruijn sequences, as they 


always are if S is a de Bruijn sequence. 


TABLE 2.5 


de BRUIJN SEQUENCE OPERATORS 








n= 3 (Canonical Form) 

S = 00011101 = HOODIES 1 

rS = 10111000 = 00010111 

S = 11100010 = OOOO 

rS = 01000111 = 00011101 
n= 4 (Canonical Form) 
S = 0000111101001011 = 0000111101001011 
rS = 1101001011110000 = 0000110100101111 
S$ = 1111000010110100 = 0000101101001111 
rS = 0010110100001111 = 0000111100101101 


If as in the case Wer m= 5 im Table 2.5ea sequence 
S 1S equivalent to its reverse complement, then S is termed 
a RC sequence. This applies to both de Bruijn and non- 
de Bruijn sequences. The following theorem appears in a 
paper by Etzion and Lempel [Ref. 6]. 
Theorem 1: A sequence S is a RC sequence if & (S) is 


even and S = {X, rX} for some X where X is a binary string. 


(hes 





F, LINEAR RECURSION FORMULA 
we brit j|n sequences ssenerated by a nonlinear function 


can also be generated by a linear function f (x, arereeNae), =) xX 


n O 


having the linear recursion formula 


: stage of the 


Were each te determines a tap to the j° 
register. The linear generator polynomial must be 
ie + ine for some positive integer k, due to the periods of 
de Bruijn sequences C2) and results from theorems given 
by Golomb [Ref. 7: pp. 27-43]. 

1. Pascal Triangle (Mod 2) 

The values for each a; are summarized in the Pascal 

triangle (mod 2) Table 2.6 for various values of n. They 


are the coefficients generated by the binomial expansion of 


ix + i)* MOG... 


TABLE 2.6 
PASCAL TRIANGLE (MOD 2) 


k = 0 1 
1 i 
2 1 0 1 
3 i i ert 
4 iL) en 
: ie 1 (x+1)° 





2. Pascal iriangle EGoperties 


(a) The a. sequence is palindromic, that is 


J 
= 4, = 1 ies a, = 4,_y> 42 = 4y-2> Cue. 
(b) The number of nonzero coefficients is even, 
n 
mat as + a. = 0 (mod 2). 
jo ? 


id 





Pt eer ee ior de BRUIJN SEQUENCES 


Pee DEPINITION OF COMPLEXITY 

The term complexity indicates a general measure of the 
predictability of a sequence. Various definitions of com- 
plexity are in use, but this paper will use the one given 
by Chan and Games [Ref. 1]. Their definition of complexity 
is the length of the shortest LFSR required to generate a 
Sequence. They apply their results primarily to de Bruijn 
sequences. This definition is also the one used by Herlestom 
[Ref. 8] to look at two different shift register cycle es 
erators. For a de Bruijn sequence S, let C(S) denote the 
Complexity of S. ° 


Previous results by Chan and Games [Ref. 1] establish 


n-l n-l 


lower and upper bounds of 2 +n and 2 respectively for 

the complexity of a de Bruijn sequence. Thus, the length 

of a LFSR is almost as long the sequence as opposed to a 
n-stage NFSR which generates the sequence. Table 3.1 serves 
fmumaeceddy sreterence tor the upper and lower limits of 
Complexity for 3<n<8. The upper bound is known to be attained 


for all n and the lower bound for all n<6. 


B. COMPLEXITY ALGORITHM 
A fast method for determining the complexity of a de 


Bruijn sequence was developed by Games and Chan [Ref. 9] 
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Txoie . 1 


bite is OF COMPLEXITY 


Span Lower Bound Upper Bound 
(n) (20 *+n) = (27-1) 
3 7 7 

4 12 15 

5 2 Sul 

6 38 a8) 

i gs ray 

8 136 255 


mers represenecea by the flow chart in Fig. 3.1. Let S be 
a de Bruijn sequence then 2(S) = 2". Further let S = A:B = D 
where A = fa, a, **> ax’ and B = {b, bs «°° Daal’ then 
Pee ta, a, °** a en-1 1 Po °°? Pon-1?- To serve as a guide 


maeexample for n = 5 1s now presented. Let the de Bruijn 


sequence S = {11111000001000110101100101001110} having length 


pueae 2 
A 1111100000100011 £ = 16 
B 0101100101001110 CaO 
A+B 1010000101101101 Cee 10 
A 10100001 es. 
B O10 dal,0 1. C = 16 
A+B 11001100 C = 24 


“9 





A 1100 Q = 4 Since A@B = 0 





B 1100 C = 24 
A+B 0000 C = 24 C is unchanged 
A al y= v2 
B 00 C = 24 
A+B 11 C = 26 
A 1 oe Add I since & = 1 
B i C = 26 G§ A= 1 
A+B 0 ee = 26 + 1 = 27 


C. COMPLEXITY DISTRIBUTION 

iitesmiature Orethe distribution of complexity of de Bruijn 
Sequences is the primary interest of this paper. Let 
a (c,n) denote the number of de Bruijn sequences having com- 
Deexity C and Span n. The complexity distribution of de 
Pan Sequences of span n for 3<n<6 are listed in Table 
See Unt ortunmacely the large number of de Bruijn sequences 
gor n>/ C25 does not allow an exhaustive examination of 
eicmr COMpPLEXITLe’s . 

Examination of this data led Chan and Games to conjecture 


that for n>3 a (c,n) is congruent to 0 modulo 4. 


Dew Pee eClS OFPOPERATORS ON COMPLEXITY 
Gin eteal. (Ret. 1) pnoved that a{c,n) for n>3 was con- 


sruent to 0 modulo 2. By showing that for each de Bruijn 
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No Yes 


Flowchart for Complexity Algorithm. 
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ABLE SS. 2 
COMPLEXITY DISTRIBUTION 


e Gale 5} 


2 


aac 50) 


8704 
18096 
34224 
67700 

L209 2 
IMSS) EIA, 
3 Oras 
1041252 
2090716 
ANG 2 552 
8342176 


16692832 
557 oi Z0 


on Cee 


64 
180 
224 
448 

1024 





sequence having complexity c, there exists another different 
de Bruijn sequence S also having complexity c. Consequently 
they showed that the complement operator does not change the 
Eomprexrty of a de Bruijn sequence. Using a different 
approach, the next theorem will show that the complexity of 
a de Bruijn sequence is also the complexity of the de Bruijn 
sequences defined by the reverse operator and by the reverse 
complement operator as well. 

Lemma 3. For S = {s, s, *** s ._,}, a de Bruijn sequence, 
mrs) = c iff E a. S, = 0 for every S5 where each a, comes 
from the oth = of the Pascal triangle (mod 2) and c is the 
smallest integer such that the above fen Oilis,. 


Proof. S is de Bruijn with C(S) = c iff 2S fF(s.i4° 


S.,.) for every s; e¢ S, which is equivalent to s, = 
c 


> a. s.,. from results by Golomb [Ref. 5: pp. 27-43] 


+ J 173 

j=1 Cc 

Wemcaoned carlier. In turn s. = 2 a; s.,. = 0 from proper- 
1 520 py 

ties of the Pascal triangle given in Section II. Q.8 -D. 


This lemma is used in the following theorem to show that 
complexity 1S preserved by various operators on de Bruijn 
sequences. 

Theorem 2. If S is a de Bruijn sequence having complexity 
C, then C(S) = C(S) = C(rS) = C(rS) = c. 

Pyoot. @ ort I (reverse opérator ) 

Let S. S; 26° S. be an arbitrary (c+l) long segment 


1+] 1+Cc 


of S which satisfies 


OS 








e 
i Y Y Sa: 
(1) a A, Sing = 0 


aekemna 5. Since the = sequence is palindromic as shown 


ler , cee Ss. S. a (c+l) long segme fe 1sfi 
earlier, ee i41 54 ( ) g gment of rS satisfies 


fomacion (1). Thus C (rS) = c. 


Part II (complement operator) substituting S; 5S ie 
> 


the c+l long segment of S, into equation (1), the results are 


S a 


es. + a. Ss. + ee + 9 ~ Ss 
1 1+c-l Cc ~1+c = 


c-l 


1+] a) iia Ge=l Se +1) mea 


(s. 


@ e 7 
i+c +1) = = a. S&.,- + fF a- (mod 2). Since there are 


0 J 
e 
an even number of nonzero a;S, x Be = 0 and equation (1) 
j=0 
is satisfied. Therefore, by Lemma 3 C(S) 


om 

Part III (reverse complement operator), by Parts I and 
II, it follows immediately that C (rS) = c. Thus C (S) = 
fens) - C (Ss) = C (rS) 


Q.E.D. 
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ier “PSO TH TABLES 


Linear and nonlinear feedback functions were used in 
conjunction with n-stage shift registers to generate de Bruijn 
Beaquences in Section Ii. The input Xg was generated as some 
function f of the current state X, Xo °°" X, that 1s 

(1) Xg = if (x) X, °°" x 
or equivalently writing f as in (2) when specializing for 


cycles only functions. 


t2) 2c > ee (x cee Xp: 
Extensive use is made of equation (2) since results for 


cycles only are desired. 


feeerULL TRUTH TABLE 

The full truth table is a useful way of listing the 
functional value Xo from equation (1) for each possible 
State of the shift register. The functional values in Table 
4.1 is an example of a full truth table for a de Bruijn 
sequence (00011101). 

Examining the full truth table it becomes apparent that 
the information in the lower half (below the dashed line) 
1s redundant, since the values of f (x, X5 Xz) fOr ez we el 
in the lower half are merely the complements of the values 
ome t (x) X5 Xz) for Xz = O in the upper half for the identical 


values of xX X>- This suggests that the upper half (or the 


Bis 


| 
| 





Current State 


TABLE 4.1 


FULL TRUTH TABLE OF (00011101) 


Functional Value 


ale 


IL 


=f (x4 X5 Xz) 


er Sewer eee 


(X3 X27 Xy) 
0 0 0 
001 
0 1 0 
011 
1 0 0 
t ‘oe 
Vik @ 
1 a 


lower half) of the truth table can display all the informa- 


tion for a function generating only cycles. 


B. HALF TRUTH TABLE 


Focusing on equation (2), since x 


1S "always | inthe 


lower half this can be viewed as complementing the functional 


values for g (x4 cee Xn-1) in the upper half. In this manner 


the half truth table in Table 4.2 is constructed for the same 


de Bruijn sequence (00011101) to aid comparison with Table 


4.1. From this point on the term truth table will always 


mean the upper half truth table with ae 0, unless noted 


otherwise. 
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TABLE 4.2 


HAP eee ABLES OF (00011101) 





The string of functional values in the truth table is 


defined to be the generator G = {g, g, ee g } where 
0 1 pe ts 

the subscripts of g are the decimal equivalent of 

ey” (CX, C- Both the weight and parity of generator G 


will be of interest as the investigation of this paper con- 
tinues. Whenever the parity or weight of the truth table 
iS mentioned in this paper, it is referring to the generator 


Of that truth table. 


C. BRO TRUTH TABLE 

The zero truth table is defined as the truth table in 
which all functional values are zero. The zero truth table 
bormm > 3 Contained in Table 4.3 is basically the output of 
a Circulating register which generates pure cycles. 

Pure c€yeies are defined as cyclic permutations of the 
Cfigindi  Staue. ihe zero truth table in Table 4.3 generates 


the four pure cycles shown in Fig. 4.1. The number of pure 


S/ 





fee 4.3 


eno TR TABLE FOR n = 3 


Xo Xp 8 (CX Xp) 


0 0 0 
Ory i 0 
1 0 0 
se 0 


cycles Z| has been shown [Ref. 7: p. 120] to be 


n 
Ze 5 6 (a2 4. 
d|n 


7 ae 
a : 

Eulers totient function ¢@ (d) is the number of fractions 
of the form */d where l<a<d and */d is in lowest terms. The 
Summation iS over all positive c*integers d which divide n, 
denoted (d{/n). It is further shown [Ref. 5] that Z- is 
the minimum possible weight truth table that can generate a 
de Bruijn sequence. 

In thewmext section, the proof of theorem 6 gives an 
example of how pure cycles can be joined in a way to create 
a de Bruijn sequence that has a minimum weight truth table. 


Table 4.5 gives a listing of Zn and Zit nen SL<ms / 


D. ONES TRUTH TABLE 
The ones truth table assigns a functional value of 1 to 
ewery entry of g (XxX) cee X,-1)- The ones truth table depicts 


the output of a complementing shift register. The input is 
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O1ll 110 


i 
| 


Paige 4.0 Pure Cycles jn=—35) . 





always the complement of the output, i.e. f (x4 30% x) = 
i+ Xs The ones truth table for n = 3 is shown in Table 
4.4. The complementing cycles from Table 4.4 are depicted 


ime Fis. 4.2. 


TABLE 4.4 
OMeESeTRUIH TABLE (n = 3) 





* 
The number of complementing cycles Z, [Ref. 7] is 
= 2d 
| ae jemeaye 2°". 


e. 
—z on 2d|{n 


* 
Since summation for Z is over all 2d (even) numbers which 


* 
divide n, for n odd Zo = an 


2 


Accordingly Fredricksen [Ref. 5] shows truth tables of max 
weight that can generate a de Bruijn sequence is 
gn-1 - 


= ap Le 
n 


= * 
ive values. of cae and 2” be Zo + teare given in Table 4.5 


for i<n<7. Also included are the number of de Bruijn se- 
quences having truth tables of maximum or minimum weight 


from a listing by Fredricksen [Ref. 5]. 
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TAB Eee 5 
PURE CYCLE COMPILATION 





# de Bruijn # de Bruijn 
sequences Sequences 
with Min er with Max 

Wt Truth Table z.* 2 -Z*4+1 Wt Truth Table 
1 i i it 
it 1 2 it 
2 Z 3 2 
Ie 2 7 3 
De Se 4 13 a 
2 ets 6 27 nol 
Paes -5 13. 10 55 cos 
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Fig. 4.2 Complementing Cycles (n=3) 


Ee “erreECTS OF CHANGING THE TRUTH TABLE 
If the ones truth table in Table 4.4 is changed so that 


Ca, |) == 0Svice 1, tle ruth table of Table 4.6 is produced. 
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TABLE 4.6 


TRUTH TABLE (TABLE 4.4 MODIFIED) 


Xo Xy g (Xx; Xp) 
0 0 1 
0 1 0 
1 0 1 
il 1 


ieroechanse in the truth table effects the cycles of Fig. 
4.2 by joining them as depicted in Fig. 4.3. Indeed the 
change produces a Hamiltonian path for the de Bruijn sequence 
eee LO111). 

Suppose the truth table of Table 4.6 is now changed so 


e 


meee (00) = "0 vice 1. This produces the truth table in 


& 


apace 4.7, and 1ts cyclic composition in Fig. 4.4. 


TBE o. 7 
TRUTH TABLE (TABLE 4.6 MODIFIED) 
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000 


001 100 
010 
101 
O1l 110 
Bul 


Pies  CYCclic C@mposition of Table 4.4. 
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000 


001 100 
010 
101 

011 ; 
Ae 


Fig. 4.4 Cyclic Composition of Table 4.7. 
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Mmiesnwcianve ine tne truth table causes the cycle to split 
mmeo two cycles. Suppose (x), eee x4) 1S on element on upper 


of two cycles as shown in Fig. 4.5 and that its successor is 


Pegs. 4.5 Iwo Cycle Graph. 


Xy-y °°° «Xz «ds as 1g (xy BOG scar) =e We tae truth table is 
now changed so that f (x4 ae ie Cemaen XxX °° XxX, 0 
1s the new successor for a Retaining the cycles 


only condition f (x, *** x x,,) must now be 1. So the 


successor of ee Oy LS OW ett XY LL .Ome of 
two things can happen how depending on whether ee °° 0 


Bogon the supper or lower cycle. 


46 





pars t 1f (X41 cee x0) 1s on the upper cycle, then that 


Seles !s Split a5 Shown in Fig. 4.6 forming three cycles. 


~ 


Go XG 
Gon KGa 


ae Soom Cycle sSplitting. 


Otherwise if (X41 2 OS x10) ue on the lower cycle, then 
those two cycles are joined as shown in Fig. 4.7. 

Clearly the parity of the cycles changes for each change 
in the truth table. The following theorem from Golomb 
[Ref. 7: p. 122] shows that the parity of the number of cycles 
and the parity of the truth table generator are equal. 

Theowem 3. For n>z, the parity of the truth table ils 
equal to the parity of the number of cycles of the truth 


table. 





Since a de Bruijn sequence 1S one cycle by definition, 
the following corollary of theorem 3 is given. 
Senomanry. iltewa sequence is de Bruijn, then the parity 


Pomme tl! Of that sequence is odd. 





fgewt,/ Crete Joining. 


Be sPECIAL NUMBERS 

Analyzing the effects of de Bruijn sequence operators 
on their truth tables, certain binary numbers are found to 
have special properties which can be exploited in later 


analysis. 
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1. TypemRe Numbers 


If X = X, "°° X, 1s a binary number, such that 
eX, TX, ey, then X 1s called a RC number. 
Lemma 4. For n even there are exactly 22 type RC numbers. 
Proof. For n even, if X is an RC number, then 
X, *** X. = X_ e** x, by definition. Therefore, 
1 n n | ala 
X) = X,> X, = X,_]> etc. Thus xyor xX, can be written as 
o ee ; in ee : : 
X1X5 XyXnXoX] - Since x positions can be tilled, in 
22 


either of two ways, the total number of RC numbers for n 


Cven is yo 
O02, EmD:, 
There are no type RC numbers for n odd, since this 


would require that x = Xne] which is impossible. 


= 


Higol 
& lype R Number§ 
If X = 





xo X, 1s a binary number such that 


eek, = XTX,» then it is called a type R number. 


— 


tO 


Lemma 5. For n even there are exactly 2° type R numbers, 


2 





and for n odd there are exactly 2 type R numbers. 
Proof. The proof for n even parallels that of Lemma 4 


mem@amwillwaot be repeated. For n odd, -since ee ae eT 


ee 


2 
there are n+l positions to be filled in either of, two 
bese Tneretore, for n odd there are exactly 2 : type R 
numbers. 


Ore 
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A useful way of organizing special numbers and other 
truth table values to assist in evaluating effects on the 
generators for various operators is shown in Table 4.8 for 
n= 5 and n = 6. The numbers are the decimal equivalent of 
tme binary (n-1)-tuple, ee from the feedback function 

=o. “n= 25) a Xn-1? 
meen also corresponds to the truth table state. 

For example for n = 5, if f (2) = 1 for some sequence 
eueemen: +f (11) = 1 for r S, f (4) = 1 for r S, and f (13) = 1 
for S. In general, for some operator 8 and some number a, 
if f (a) = 1 for some sequence S, then f (fa) = 1 for B S. 
Special numbers have properties in addition to this. If a 
is a type RC number and f (a) = 1 for S, then f (a)= £ (ra) = 
1 for r S. Similarly if y is a type R number and f (y) = 1 
Pores, then £ (y) = f (ry) = 1for rs. Thus, if we want to 
shown S is a sequence such that S = r S, then if f (a) = 1 


in S, then f (r a) = 1 in S also. 


G. GENERATOR ANALYSIS 

Analysis of the generators for the de Bruijn sequences 
S, r S, S, & r S in Table 4.9 appears to indicate the gener- 
ator G of sequence S denoted G (S) is the reverse of the 
pemerator G (S$). Notice also that G (r S) = r [G (r $)] 
fer os. Althewoh Table 4.9 is a@n €xample of a singhe sequence 
Dec 1S@met aneise lated incident; im general the result 
G (S) = r [G (S)] holds and is proven in the following 


theorem. 
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TABLE 4.8 


TRUTH TABLE ANALYSIS CHART 


5 S es bee 
S S ios) rs 
cS gs S S 
£5 rs S S 
0 i 0 15 i« Truth table states 
1 14 8 fi Type R numbers 
2 ts 4 11 (0FGe 9,15) 
5 12 LZ 3 
Type RC numbers 
5 10 10 5 
(Soe 012) 
6 9 6 9 


«Truth table states 
Type R numbers 
GOR LO 51 } 
Type RC numbers 


(None) 
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TABLE 4.9 


GENERATOR ANALYSIS FOR de BRUIJN SEQUENCES 
(r S # S) 


S 11111001000111011010011000001010 


nI 


11111010100000110111000100101100 
es) 11111010100000110010110111000100 
rs 11111001101001000111011000001010 


ne | GC GSa G (S) G ae) Ge (rs) 


? 


1 1 J, I. 

0 IL i 0 

1 0 0 1 

/ | 1 0 0 il 
aye RC 0 i il 0 
0 il i 0 

C) 0 == 0 1 
Type R 0 1 i 0 


nner ay 
0§ 15 


ee oe 
_ co ~“ mi = Iwill nm KF oO 
}— 
S 
} 
SS 


0 1 0 0 1 

i 0 0 1 
lly 0 i 1 0 
15 0 i 1 0 
14 iL 0 0 1 
15 1 1 1 iL 
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Theorem 4: If S is a de Bruijn sequence with generator 
G (S), then S is a de Bruijn sequence with generator 
fic (S)]. 


Proof: Let G (S) = {gy > ime ss > & } then the full 


on-1 


meme cable for S is given in Table 4.10. 


TABLE 4.10 


FULL TRUTH TABLE FOR S 





Iie x 


n? ee ee )eiseamined) long string 


plea 
in S, then mapas’ ee Xp X, Xo °° Xx 1s a (merl]) 
long string in S. 

Therefore we can transform each line of the full truth 


table for S by taking its complement to get Table 4.11. 





TABLE 4.11 


FULL TRUTH TABLE FOR S (NON STANDARD) 


ST Sell il f (Xx); 2 x) 
E ji if ve. it gy 
1 if ip cee 0 g5 
1 0 0 vee OT UG 
on il 
0 1 1 cee 1 gy 
0 1 1 “se 0 g5 
0 0 0 seks 0 Cire. 
; on 1 


Rearranging the rows from top to bottom, the TT for S is 
placed in standard form as depicted in the Table 4.12. 

Further examination of the generators in Table 4.9 shows 

the structure of G (S) to be remarkably similar to the 
Seructure Of G (r S), with the exception that the functional 
values of type R numbers (6 §& 9) are transposed. If, however, 
these values should be the ae then clearly G (S) = G (r S) 
or S = r S. Indeed, this occurs twice for n = 35 and 64 times 
for n = 5, and will be shown to happen for all odd n>3. 

This cannot happen for n even, because there are no type 


RC numbers in the truth table. 
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TABLE 4.12 


FULL TRUTH TABLE FOR S 





Olewely G (S) as r [G (S)] 
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V. RESULTS 


Examining Table 3.2 it can be quickly determined that the 
number of de Bruijn sequences of given complexity for 3<n<6 
is 0 (mod 4). Games and Chan [Ref. 1] made this conjecture, 
but were able to show only that the distribution a (c, n) 
was 0 mod 2. Their result followed by showing that for a 
peven de Bruijn sequence S of complexity c, there existed 
a unique de Bruijn sequence S also having complexity c. 

In Section III it was shown that reverse and reverse dual 
Operators preserve the complexity of de Bruijn sequences as 
well. 

If each of these new operators produced unique de Bruijn 
sequences, then it would be easy to show that the numbers 


0 (mod 4) since for every de Bruijn sequence S&S 


omc, n) 
there would also exist S, r S and r S all de Bruijn and 
unique having complexity c. As mentioned in Section IV, the 
sequences r S and r S are not always distinct from the 
sequences S §& S respectively; however, for n even the next 
theorem shows that all four of these sequences are distinct. 
In addition it will be shown that RC sequences occur for 

all n odd, and a lower bound for the number of RC sequences 


will be established for all n and investigated for n = 7. 
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Peeeeommc, mf) = OFMODULO 4 FOR n EVEN 

Etzion and Lempel [Ref. 6] give a proof that a (c, n) 
mmcomeruent to 0 mod 4 for c n even. In the proof for n 
even given here a different approach is taken. 

iieeremss. For even n>53, the numbers a (c, n) = 0 (mod 4) 
for a de Bruijn sequence S of span n and complexity c. 

maoot. Prom theorem 1, it is known that for a de Bruijn 
sequence S of complexity c, that S, r S & r § also have 
complexity c. Games and Chan showed the uniqueness of S, 
Bommel! that rem@ins 1S to show that r S$ is distinct from 
S, or equivalently that r S is distinct from S for n even. 


Suppose on the contrary that S = r S. Then consider 


an n-tuple 0a,_, °** a, in the truth table such that 

ig (a) ei) 0) = 1. Then Oa 1 cee ay lietssa nde long 
string in S, and o ay oes Ua l1 is a n+l long string inr S. 
heeS-= r S 


then 0 ay oo oa! 1S lone string in S, 


which implies 
that 1 a, *** a,_, 0 is a n+l long string in S, and 
Pm@eretoresf (a) °° ane] 0) = 1. This must hold for every 


0 such that f (xX) © *senuexX toe = Is, «Patrane 


zak Xn-1 


these vectors, a, *** a, _, and a,j °° a 


16 asa | 
1? two at a time 


pesults in@an ewen number of 1's in the, truth table. So 
there must be a vector which pairs with itself as the parity 


Steoncsmisaodd) Bue then Xe °° Xp Xz ee Xyez TS 2 


4 


| 





contradiction since x for n even. Therefore for 


me +1 
2 y 


Mmeeven os - r 5S and r S # 5S, so r S and r § are distinct 








oma a (c, n) = 0 mod 4 
5/2 De 
B. RC SEQUENCES EXIST FOR ALL ODD n>5 
The results a(c,n) = 0 (mod 4) for all n are not obtained 


as hoped for. For n = 3 and n = 5 RC sequences are known 
Bo exist, but it is not known if this is true for all odd 
n>3 or just some. The following theorem shows that RC 
meqecences exist for all odd n>3. Recall that RC sequences 
are those sequences S for which S = r S or equivalent by 
Ge= + S since r (S) = r (r S$) yields S = r S. 

Lemma 6. The sequence formed by joining all pure 
cycles of n-tuples is de Bruijn. 

Oot newer: clemenmts Of Yeach pire cycle for a 
given n are easily composed and arranged in order by weight 
according to the algorithm in Appendix A. Since the pure 
cycles contain all n-tuples exactly once, all that remains 
is to show how each pure cycle can be interconnected, thus 
Forming a dem Bruijn sequence. For each cycle of wt a>l 
there exists a least element "2X+1" (necessarily odd, if 
it were even (2X) then "X" is on the same cycle and smaller) 
whose predecessor "X'"' has wt (a-l). Since the pure cycles 
ieee Ordemecampy weleght, it 1s clear that the least element 


Oipeacheemrescycle of wi a>) cen be joined to a cycle of 


28 





‘wt (a-1) by changing the truth table for its predecessor, 
l.e. g(X) becomes 1. In this way all the pure cycles are 
Meweinmterconnected, and the sequence is de Bruijn. 
Ore. D., 

fiteemem O- Forvall odd n>3, there exist RC de Bruijn 
Feauences. 

Proof. An RC de Bruijn sequence can always be composed 
immeene tollowing way for n>3. First construct a light cycle 
(LC) for those pure cycles of wt Foe identical to the 


— 2 
method used in Lemma 6. Note that when the cycles of 





wt — avemcmm@eamire cycles of smaller weight the truth 


table is changed at positions whose weight is less than 
a , For every X such that f (X) = 1 in LC, induce 

ghar Jt 
fe) = 1 Om the pure cycles of wt (n-a) > 2 , forming 


mmmeavy cycle (H€). Note that HC is the reverse complement 
of LC. All of the positions changed so that g (X) = 1 on 


HC are of weights bigger than or equal to a . 





Join LC and HC by using a type RC number, y, which exist 
by Lemma 4, since for n odd the truth table states have 
even (n-l) lemeth. The RC number y haS weight a ag cuca elke 
im (1) > leeiem joins LC @nd HC forming one sequence. By 
Theorem 1 the sequence is RC, and by Lemma 6 the sequence 


fede Bruin. 


Q.E.D. 


ee 





C. LOWER BOUND ON NUMBER OF RC SEQUENCES 
re tieemen 
A weak lower bound on the number of RC de Bruijn 
sequences is established. 
Theorem 7; A lower bound on the MUMbete or dew Bruzjn 


Eo) 
Pemoequences for odd n>3 is A» 2 » wheme A 1s the 


mMumper Of possible different interconnections for pure 


n-l 
Sy 





cycles of weight dove teuminea by the Best method 


{Ref. 5]. 

Pasoot:> Restricting attention to those cycles of 
wa 
motemethod gives the number A of possible different inter- 


(since the other connections are induced), the 





connections for, the Dlme Gree s) (OL weight<35* -n-pince 


g 


Ehere exist 2 type RC numbers, this gives 2 : ways 


of Joining Pedtemerent Pairs of cye€les for a total of 
A+ 24 possibilities. 


Q.E.D. 


¢. Example for n = 7 


iigememme cycles Otewt<3 sand the adjacency matrix 
Com wabeweavenminelables 5.3 and 5.4 respectively. The 
value of the determinant of the adjacency matrix is 38,880 
for any diagonal element. This gives a lower bound of 


38,880 X 8 = 311,040 for the number of de Bruijn RC sequences. 
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-_ 7 


irate 5.1 


PUREMCYE@ES TOR N=7 G wt < 3 


+ 8 16 oy 
2 24 48 26 
20 40 80 55 
28 26 112 a7 
36 WeZ Ly 54 
44 88 49 98 
S14 104 81 55 
76 Ue 20 100 
84 41 82 S)/ 
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Hebi Os. 2 


ADJACENCY MATRIX FOR n=7 
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ver, CeNdJECTURES 


Although a complete proof for a (c, n) = 0 (mod 4) 
has yet to be given, a technique for grouping RC sequences 
in a 0 (mod 4) fashion will be demonstrated for n= 5. 
mewevyer, the technique fails to preserve complexity in 
every case. However, for n = 5 the complexitites produced 
were congruent to 0 (mod 2) and when the complement of 
each sequence is added a 0 (mod 4) distribution results. 
It remains, however, that the complexity of the generated 
Sequences are not determined apriori. The technique pre- 
sented is interesting in its own right, and it is hoped 
that the interested reader may be able to apply it towards 


@eeorution of the 0 (mod 4) distribution. 


A. RC SEQUENCE GENERATOR TECHNIQUE 

Table 6.1 shows a listing of four RC sequences and their 
generators. The technique used was to choose in turn each 
of the type RC numbers to have a functional value of l. 

Note that the complexity of Sz and Sc are both 23, 
and §S 


while the complexity of S are both 29. At present 


10 lez 
there is no way of Knowing which pair of RC sequences will 

have the same complexity, except that it appears C(seee(s 
iffall four RC sequences have the same complexity. A listing 


of the generators for each of the 32 pairs of RC sequences 
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Gi GW Wh 


C2) 


TABLE 6.1 
RC GENERATOR TECHNIQUE 
C(Sz)=23 C(S-)=25 C(Sj9)=29 


eS.) G(S.) Sai) 
r 1 1 
1 1 1 
0 0 0 
ilies 0 0 
1 1 1 
0 a 0 
1 1 1 
1 1 l 
1 1 1 
1 1 1 
0 0 1 
0 0 0 
0 0 0 
1 1 1 
1 1 1 
1 1 1 


indicates RC number. 

11111001010001001101110101100000 
11111001011000001101110101000100 
11111001010110000011011101000100 


= 11111000001101110101100101000100 
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GOS 
G(S 
A: 


1 


122729 


12? 





is prepared in Table 6.2. These are listed in order by 

W(G) and then by complexity as far as the groupings allow. 
Keep in mind, that the complement of each RC £eGBice 1s not 
included in the table. A close analysis of this table then 
Shows a 0 (mod 4) grouping for the number of RC sequences 


of complexity c and span n = 5S. 


B. PROBLEMS 

For n = 7, there exist 6 type R numbers (other than 0 
& 64) and 8 type RC numbers in the truth table. This large 
number of possibilities will create considerable problems 
in extending this idea to n = 7 and beyond. 

Using this approach a large scale computer based analysis 
is essential to any further investigation of RC sequences 
Beem 7. Considering the number of de Bruijn sequences for 


7 7 


n= 7 to be 2 or 1.42 x 10 an algorithm for generating 


only the RC sequences is crucial, since generating ror 
de Bruijn sequences is not technically feasible. Though 
Theorem 7 establishes the lower bound of 311,040 RC sequences 
for n = 7, this is a very weak bound and a presumably still 


6 


conservative estimate would be on the order of 10° actual 


RC sequences for n = 7. 
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100 


100 


101 


100 


101 


100 


iglich 


OL 


COfF CCF CF oF oo Oo O}F 


or 


TABLE 6.2 


RC de BRUIJN SEQUENCE 0 (Mod 4) 
GROUPING FOR n=5 


(S) WEG(S) ] C(S) 
1001 Q 101 7 ZS 
IAS 
"| 
27 
ey, 
ae) 
a 
78 
aL 
oe. 
SL 
oF 
cat 
SM 
31 
al 
og! 
51 
pl 
cal 
31 
Sal 
a 
Sl 
Q1l Bl ns 
2 

dil. (es, 

lig LS. 
a 11 25 

lag A 
ZT 
rae 


OOF OFDM 
Ol © 
j- 


LOmet 


oj © 


o|e 


0010 


— 
Cj rk 
© 
= 
_ 


COOlr 


0010 


Se) ey 
oH 


16P1 


oo 
OfH Oo 


OlH 


1011 


O|re 


Ole © 
WOMOWMWWW WHO NYA YNYNNYN YN YIN 


|r 


Poe 


Ol} oO 
-_— 
_— 


POL be OL 


Ol HH 


jr 
Ole © 
br 
4 
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APPENDIX A 
ALGORITHM FOR GENERATING LEAST ELEMENT 
Pereeeeery PURE CYCLE 
This algorithm is designed to produce the least element 
of each pure cycle. The least elements are initially in 
@eereasing order. A simple rearrangement will group the 
least elements by weight in decreasing order if necessary. 
mie algorithm to be presented is an adaptation of the 
8-operation discussed in a paper by Fredricksen and Maiorana 
[Ref. 10], to generate a lexicographic list of necklaces. 
ALGORITHM: Begin the 2 operation with the zero n-tuple,. 
feet 1S the first least element for the zero pure cycle. 
2 Operation: 2 (x) coe x = (Y) cee ee, = Y 
1. Find the largest subscript j such that 
x = 0 and xX, = Loree). 
2. Form x, xX, °° Xi o4 1 where Xi = 1 
3. Repeat x); xX, °° X51 until n numbers are produced. 
iwi = £) “<enusmme Ineeger t then let 
YexyooexX; ylxyeeexy yleeexgee xy yl 


11) If n>tj then finish with a a 
J 
then let Wa eee 9 5-2 % °° ** 


4. Y is a least element iff n=tj 


-t.° 
ais, 


5. Repeat beginning at step 1 until the final least 


element ca is reached. 
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TABLE A.1 


2 ALGORITHM FOR N = 7 


Wt 
0000000 0 
0000001 1 
0000010 * 

0000011 2 
0000100 * 

0000101 2 
0000110 * 

0000111 3 
0001000 * 

0001001 ae 
0001010 * 

0001011 3 
0001100 * 

0001101 3 
0001110 * 

0001111 4 
0010010 * 

0010011 3 
0010100 * 

0010101 3 
0010110 * 

0010111 4 
0011001 * 

0011010 * 

0011011 4 
0011100 * 

0011101 4 
0011110 * 

0011111 5 
0101010 * 
0101011 4 
0101101 * 

0101110 * 

0101111 5 
0110110 * 

0110111 : 
0111011 * 

0111101 * 

0111110 * 

il alt 6 
Liege 7 


* Unacceptable 


68 





tavhewn.) illustrates the algorithm for n=7. 


Arranging 


mmeeader Dy Welght, the least element for each pwre cycle is 


tabulated in Table A.2. 


PEURESOVGEE 


coe 


“ION 


TABLE A.2 


LEAST ELEMENTS (n=7) 


0000000 
0000001 
0000011 
0000101 
0001001 
0000111 
0001011 
0001101 
0010011 
0010101 


Creat 1 
OC u 11 
Cree 
0011101 
0 e021 
0011111 
@i01111 
OnenOHe 1 
OF isi 
Moll 2 1 
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